Dynamical traps and chaotic advection in a meandering jet flow 
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We continue our study of chaotic mixing and transport of passive particles in a simple model of 
a meandering jet flow [Prants, et al, Chaos 16, 033117 (2006)]. In the present paper we study and 
explain phenomenologically a connection between dynamical, topological, and statistical properties 
of chaotic mixing and transport in the model flow in terms of dynamical traps, singular zones in the 
phase space where particles may spend arbitrary long but flnite time [Zaslavsky, Phys. D 168—169, 
292 (2002)]. The transport of passive particles is described in terms of lengths and durations of zonal 
flights which are events between two successive changes of sign of zonal velocity. Some peculiarities 
of the respective probability density functions for short flights are proven to be caused by the so- 
called rotational-islands traps connected with the boundaries of resonant islands (including those 
of the vortex cores) fllled with the particles moving in the same frame. Whereas, the statistics of 
long flights can be explained by the influence of the so-called ballistic-islands traps fllled with the 
particles moving from a frame to frame. 
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I. INTRODUCTION 

Methods of the theory of dynamical systems are ac- 
tively used to describe advection of water (air) masses 
and their properties in the ocean and atmosphere 

Q, Q, 3, B Q, Q, 0, Q, H, Q, Q, H, 

The geophysical jet currents, like the Gulf Stream, the 
Kuroshio, the Antarctic circumpolar current, and oth- 
ers in the ocean and the polar night Antarctic jet in the 
atmosphere, are robust structures whose form typically 
changes in space and time in a meander-like way. If ad- 
vected particles rapidly adjust their own velocity to that 
of background flow and do not affect the flow properties, 
then the particles are called passive (scalars, tracers, or 
Lagrangian particles) and their equation of motion is very 
simple 



dr 
'dt 



= V(r,i) 



(1) 



where r = {x,y,z) and v = {u,v,w) are the position 
vector and the particle velocity at a point {x,y,z). If 



the corresponding Eulerian velocity fleld is supposed to 
be regular, the vector Eq. (U) in nontrivial cases is a 
set of three nonlinear deterministic differential equations 
whose phase space is a physical space of advected par- 
ticles. It is well known from dynamical systems theory 
that solutions of this kind of equations can be chaotic 
in the sense of exponential sensitivity to small variations 
in initial conditions and/or control parameters. As to 
advection equations, it was Arnold ^ who firstly sug- 
gested chaos in the field lines (and, therefore, in trajec- 
tories) for a special class of three-dimensional stationary 
fiows (so-called ABC fiows), and this suggestion has been 
confirmed numerically by Henon 0]. In the approxima- 
tion of incompressible planar fiows, the velocity conipo- 
nents can be expressed in terms of a stream function [lol ]: 
u = —d'^/dy and v — d"^ /dx. The equations of motion 
(HJ) in a two-dimensional incompressible fiow have now a 
Hamiltonian form with the streamfunction 4" playing the 
role of a Hamiltonian and the coordinates {x, y) of a par- 
ticle playing the role of canonically conjugated variables. 
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chaotic in a number of theoretical UJ 
n, Q, 0, Q, 01 and laboratory Q, H 




Advection of passive particles has been shown to be 

lyy 

2l| models 

of geophysical jet currents. In our recent paper jl| (Pa- 
per I) we have studied mixing, transport, and chaotic 
advection in a simple kinematic model of a meandering 
two-dimensional jet flow with a Bickley zonal velocity 
profile u ~ sech^ y being motivated by the problem of 
lateral mixing of water masses (together with salinity, 
heat, nutrients, pollutants, and other passive scalars) in 
the western boundary currents in the ocean. We derived 
advection equations in the frame moving with the phase 
velocity of a running wave imposed on the Bickley jet, 
found the stationary points, conditions of their stability, 
and all the possible bifurcations of these equations which 
were shown to be autonomous in the co-moving frame. 
Under a periodic perturbation of the wave amplitude, the 
phase plane of the chosen model flow has been shown to 
consist of a central eastward jet, peripheral westward cur- 
rents, and chains of northern and southern circulations 
(vortex cores) immersed in a chaotic sea which, in turn, 
contains islands of regular motion. Statistical properties 
of chaotic transport of advected particles have been char- 
acterized in terms of particle's zonal flights (any event be- 
tween two successive changes of the sign of the particle's 
zonal velocity u). Probability density functions (PDFs) 
of durations and lengths of flights, computed with a num- 
ber of very long chaotic trajectories, have been found to 
be complicated functions with local maxima and frag- 
ments with exponential and power-law decays. 

The aim of this paper is to study a phenomenological 
connection between dynamical, topological, and statis- 
tical properties of chaotic mixing and transport in the 
meandering jet flow considered in Paper I and to explain 
transport properties by phenomenon of the so-called dy- 
namical traps. Following to Zaslavsky the dynamical 
trap is a domain in the phase space of a Hamiltonian sys- 
tem where a particle (or, its trajectory) can spend arbi- 
trary long flnite time, performing almost regular motion. 



despite the fact that the full trajectory is chaotic in any 
appropriate sense. In fact, it is the definition of a quasi- 
trap. Absolute traps, where particles could spend an infi- 
nite time, are possible in Hamiltonian systems only with 
a zero measure set. The dynamical traps are due to a 
stickiness of trajectories to some singular domains in the 
phase space, largely, to the boundaries of resonant islands 
and cantori. There are no classification and description 
of the dynamical traps. Zaslavsky described two types 
of dynamical traps in Hamiltonian systems: hierarchical- 
islands traps around chains of resonant islands Q,0,Qi 
and stochastic-layer traps which are stochastic jets inside 
a stochastic sea where trajectories can spend a very long 
time 0, 0, It is expected that classification and de- 
scription of the most typical dynamical traps would help 
us to construct kinetic equations which will be able to de- 
scribe transport properties of chaotic systems including 
anomalous ones [allll. Il2|. 



This paper is organized as follows. We start with ad- 
vection equations derived in Paper I in the frame of refer- 
ence moving with the phase velocity of a meander whose 
amplitude changes in time periodically. We compute in 
Sec. II PDFs for lengths Xf and durations Tf of zonal 
fiights for a number of chaotic trajectories and show an- 
alytically that all the fiights start and finish only inside 
a strip confined by two curves whose form is defined by 
the condition u — 0. The prominent peaks in statistics 
of short fiights (|a;/| < 27r) are proved to be caused by 
stickiness of trajectories to the boundaries of rotational 
resonant islands filled with regular particles rotating in 
the same frame (Sec. III). We call this type of dynami- 
cal traps as a rotational-island trap (RIT). Another type 
of islands, ballistic islands (filled with regular particles 
moving from frame to frame), is proved to contribute to 
the statistics of long fiights (|x/| ^ 27r) in Sec. IV. Both 
the ballistic-islands trap (BITs) and RITs belong to the 
class of hierarchical-island traps by the Zaslavsky's clas- 
sification. 
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II. DESCRIPTION OF CHAOTIC TRANSPORT 
IN TERMS OF FLIGHTS 

A. Basic features of the flow 

We take the following specific stream function as a 
kinematic model of a meandering jet flow in the labo- 
ratory frame of reference: 



= -^-otanh 



y' — a cos k{x' — ct') 
\^l + k^a^ ain^ k{x' - ct') ^ 



(2) 



where the width of the jet is A. Meandering is provided by 
a running wave with the amplitude a, the wave-number fc, 
and the phase velocity c. The normalized streamfunction 
in the frame moving with the phase velocity is 



= - tanh 



y ^ A cos X 

sin^ X , 



Cy, 



(3) 



where x — k{x' — ct') and y — ky' are new scaled co- 
ordinates, and A — ak, L — Xk, and C = c/^'q/c are 
the control parameters. Equations, governing advection 
of passive particles |(T]) in the co-moving frame, are the 
following: 

1 



C, 



+ A^ sin^ x cosh 9 
^sina;(l + — Ay cosx) 

y ~ A cos X 



(4) 



+ sin^ X 

where dot denotes differentiation with respect to dimen- 
sionless time t = ^'^k'^t' . 

In Paper I (for more details see [l^) we have found and 
analyzed all the stationary points, their stability, and the 
bifurcations of the equations of motion ^ . Being moti- 
vated by the problem of mixing and transport of water 
masses and their properties in oceanic western bound- 
ary currents like the Gulf Stream and the Kuroshio, we 
chose the phase portrait shown in Fig. [l^ among all 
the possible flow regimes. Passive particles can move 
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FIG. 1: (a) Stationary streamfunction of a meandering jet 
in the co- moving frame ([S}. The flow is divided into three 
different regimes: circulations (C), jet (J), and peripheral 
currents (P). (b) Poincare section of the perturbed mean- 
dering jet in the co-moving frame. The parameters of the 
steady flow are: the jet's width L = 0.628, the meander's am- 
plitude Ao = 0.785 and its phase velocity C = 0.1168. The 
perturbation amplitude and frequency are: e = 0.0785 and 
w = 0.2536. 

along stationary (in the co-moving frame) streamlines 
in a different manner. They can move to the east in 
the jet (J) and to the west in northern and southern 
(with respect to the jet) peripheral currents (P). There 
are also particles rotating in the northern and south- 
ern circulation cells (C) in a periodic way. The north- 

(n) 

ern separatrix connects the saddle points at Xs = 27rn 
and i/s"'' = L Arcosh ^JxJlC -I- A and the southern one 

is] 

connects the saddle points at xl — {2n -I- l)7r and 
yi"^ = -L Arcosh ^Jl/LC + A, where n = 0, ±1, . . . . 

As a perturbation, we took in Paper I the simple peri- 
odic modulation of the meander's amplitude 

A = Aq + e cos ujt. (5) 

Under the perturbation, there arise resonances between 
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the perturbation frequency w and the frequencies / of the 
particle's rotation in the circulations C. A frequency map 
fi^Oi yo), computed in Paper I (see Fig. 2 in that paper), 
shows the values of / for particles with initial positions 
{xo,yo) in the unperturbed flow. With a given value of 
the perturbation frequency and fixed values of the other 
control parameters, specific regions of the unperturbed 
fiow are destroyed. In Paper I we fixed the scaled val- 
ues of the parameters of the unperturbed fiow, the jet's 
width L = 0.628, the meander's amplitude — 0.785 
and its phase velocity C = 0.1168, that are in the range 
of the realistic values for the Gulf Stream 
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The 



perturbation frequency w — 0.2536 chosen in Paper I is 
close to the values of the rotation frequency / of the par- 
ticles circulating in the inner core of the regions C (see 
Fig. 2 in Paper I) . In Fig. [1] we show the Poincare sec- 
tion (for a large number of trajectories) of the meander- 
ing jet whose amplitude is modulated with the frequency 
w = 0.2536 and the strength e = 0.0785. The vortex 
cores survive under this perturbation, the stochastic lay- 
ers appear along the unperturbed separatrix, and a cen- 
tral jet J is a barrier to transport of particles across the 
jet. 

The equations of motion (0]) with the perturbation ^ 
are symmetric under the following transformations: (1) 
t^t, x— s-TT + x, y— !• —y and (2) t — i, x 
—X, y ^ y. It implies that the meridional transport 
(north-south and south-north) is symmetric but the zonal 
transport (west-east and east-west) is symmetric under 
a time reversal. Due to these symmetries motion can be 
considered on the cyHnder with < a; < 27r and y > 0. 
The part of the phase space with 2?™ < a; < 2n{n + 1), 
71 = 0, ±1, . . . , is called a frame. 

It is convenient to characterize chaotic mixing and 
transport in terms of zonal fiights. A zonal fiight is a 
motion of a particle between two successive changes of 
signs of its zonal velocity, i. e. the motion between two 
successive events x — u = 0. Particles (and correspond- 
ing trajectories) in chaotic jet fiows can be classified in 



terms of the lengths of fiights Xf as follows. The trajec- 
tories with \xf\ < 2tt correspond to the particles moving 
in the same frame or in neighbour frames. There are par- 
ticles moving chaotically forever in the same frame but 
they are of a zero measure. Among the particles with 
inter-frame motion, there are regular and chaotic bal- 
listic ones. Regular baUistic trajectories can be defined 
as those which cannot have two fiights with |a:/| > 27r 
in succession. They correspond to particles moving in 
regular regions of the phase space persisting under the 
perturbation, (eastward motion in the jet and western 
motion in the peripheral current) and those moving in 
the stochastic layer (trajectories belonging to baUistic 
islands). Typical chaotic trajectories have complicated 
distributions over the lengths and durations of fiights. 

B. Turning points 

As in Paper I, we will characterize statistical proper- 
ties of chaotic transport by probability density functions 
(PDFs) of lengths of fiights P{xf) and durations of fiights 
P(Tf) for a number of very long chaotic trajectory. Both 
regular and chaotic particles may change many times the 
sign of their zonal velocity x = u. From the condition 
i = in Eq. (0]), it is easy to find the equations for the 
curves which are loci of turning points 



Y±{x,A) = ±iVl + sin^c 



.XX 



X Arsech + siv?' X + Acosx. (6) 

We consider the northern curve, i. e., Eq. ^ with the 
positive sign. Taking into account that the perturbation 
has the form ([5]) , we realize that all the northern turning 
points are inside a strip confined by two curves of the 
form ([6]) with A — Aq±£. Let us analyze the derivative 
over the varying parameter A 

dY 



dA 



ACL^ sin^ X 
2D 



2 Arsech - 



1 



I, (7) 
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where D = LC^/l + A^ sir? X. If the derivative at a 
fixed value of x does not change its sign on the interval 
Aq ~ e < A < Aq + e, then Y varies from Y{x, Aq ~ 
e) to Y{x,At^ + e), and for each value of y we have a 
single value of the perturbation parameter A. However, 
there may exist such values of x for which the equation 
dY/dA — has a solution on the interval mentioned 
above. In this case one may have more than one values 
of A for a single value of y. Thus, the width of the strip, 
containing turning points, is defined by the values of Y at 
the extremum points and at the end points of the interval 
of the values of A. In Fig. [2] we show the turning points 
of a single chaotic trajectory on the cylinder < x < 2n 
confined between two corresponding curves. 
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FIG. 2: Turning points of a single chaotic trajectory on the 
cylinder < x < 2n are in a strip confined by two curves ((G} 
with A = Ao±e 

To study chaotic transport we have carried out nu- 
merical experiments with tracers initially placed in the 
stochastic layer. It was found that statistical properties 
of chaotic transport practically do not depend on the 
number of tracers provided that the corresponding tra- 
jectories are sufficiently long (t ~ 10^). The PDFs for 
the lengths x / and durations T/ of fiights for five tracers 
with the computation time t = 5 ■ 10* for each tracer 
are shown in Figs. [SK and b, respectively, both for the 
eastward (e) and westward (w) motion. Both P{xf) and 
P{Tf) are complicated functions with local extrema de- 
caying in a different manner for different ranges ofxf and 
Tf . The main aim of our study of chaotic transport is to 



figure out the basic peculiarities of the statistics and at- 
tribute them to specific zones in the phase space, namely, 
to dynamical traps strongly infiuencing the transport. 
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FIG. 3: Probability density functions of (a) lengths Xf and 
(b) durations Tf of the westward (w) and eastward (e) flights. 
The PDFs PwiTf) and Pe{Tf) are normalized to the number 
of westward and eastward flights, respectively. 



III. ROTATIONAL-ISLANDS TRAPS 

It is well known, that in nonlinear Hamiltonian sys- 
tems a complicated structure of the phase space with is- 
lands, stochastic layers, and chains of islands, immersed 
in a stochastic sea, arises under a perturbation due to 
a variety of nonlinear resonances and their overlapping 
[2^ . The motion is quasiperiodic and stable in the is- 
lands. The boundaries of the islands are absolute barriers 
to transport: particles can not go through them neither 
from inside nor from outside. In typical Hamiltonian sys- 
tems there are also partial barriers called cantori (for a 
review see j^l)- The cantori are invariant Cantor sets in 
which the motion is quasiperiodic and unstable. They re- 
semble invariant tori but have an infinite number of gaps 
through which irregular particles can penetrate. How- 
ever, the transport through the gaps is difficult since the 
total length of the gaps is finite and, therefore, they are 
very small. The cantori present strong partial barriers to 
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transport. 

Invariant curves of the unperturbed system (Fig. [T^) 
are destroyed under the perturbation |[5]). As the pertur- 
bation strength e increases, a closed invariant curve with 
frequency / is destroyed at some critical value of e. If 
the f /w is a rational number, the corresponding curve is 
replaces by an island chain, while the curves with irra- 
tional frequencies are replaced by cantori. There are un- 
countably many cantori forming a complicated hierarchy. 
Numerical experiments with a variety of Hamiltonian sys- 
tems with different number of degrees of freedom provide 
an evidence for the presence of strong partial barriers to 
transport around the island's boundaries (for review, see 
I4I) which manifest themselves on Poincare sections as 
domains with increased density of points. 

In Paper I we have found that with chosen values of the 
control parameters there exist in each frame a vortex core 
(which is an island of the primary resonance w — f) im- 
mersed into a stochastic sea, where there are six islands 
of a secondary resonance emerged from three islands of 
the primary resonance 3/ = 2w (see Fig. 3 in Paper I). 
Chains of smaller size islands are present around the vor- 
tex core and the secondary-resonance islands. Particles 
belonging to all of these islands (including the vortex 
core) rotate in the same frame performing short flights 
with the lengths < 27r. So we will call them rota- 
tional islands and distinguish from the so-called ballistic 
islands to be considered below. 

Stickiness of particles to boundaries of the rotational 
islands has been demonstrated in Paper I. It means that 
real fluid particles can be trapped for a long time in a 
singular zone nearby the borders of the rotational islands 
which we will call rotational-islands traps (RITs). To il- 
lustrate the effect of the RITs we demonstrate in Figs. [4] 
and [5] the Poincare sections of a chaotic trajectory in the 
frame < a: < 27r sticking to the vortex core and to the 
secondary-resonance islands, respectively. The contour of 
the vortex core is shown in Fig. [5] by the thick line. The 
small points are tracks of the particle's position at the 



moments of time i„ = 27rn/w (where n = 1,2,...) and 
the thin curves are fragments of the corresponding trajec- 
tory on the phase plane. Increased density of points indi- 
cates the presence of dynamical traps near the boundaries 
of the rotational islands. Contribution of the vortex-core 
RIT (Fig. |4]) to chaotic transport is expected to be much 
more significant than the one of the RITs of the other 
islands (Fig. [5]). 




3.6 



FIG. 4: The vortex-core trap. Poincare section of a chaotic 
trajectory in the frame < a; < 27r with a fragment of a 
trajectory. 




FIG. 5: The secondary resonance islands trap. A fragment of 
a chaotic trajectory sticking to the islands is shown. 

It is reasonable to suppose that RITs contribute to the 
statistics of short fiights. By short fiights we mean the 
flights with the length shorter than 27r. In Fig. [6] we 
show the part of the full PDF P{Tf) (Fig. [Sb) for the 
eastward (e) and westward (w) short flights separately. 
There are a comparatively small number of the eastward 
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FIG. 6: The PDFs for the eastward (e) and westward (w) 
flights with the length shorter than 2tv. The PDFs Pw{Tf) 
and PsiTf) are normalized to the number of westward and 
eastward flights, respectively. 
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FIG. 7: The vortex-core trap PDFs of durations of the 
eastward (e) and westward (w) flights Tj for (a) a regular 
quasiperiodic trajectory inside the vortex core and close to its 
boundary and (b) a chaotic trajectory sticking to the bound- 
ary of the vortex core from the outside. 

flights with Tf < 11. Let us note the prominent peak 
of the corresponding PDF at ~ 11 fohowed by an 
exponential decay. As to the westward short flights, there 
are two smah local peaks around ~ 17 and 21. 

To estimate the contribution of the vortex-core RIT to 
the statistics of short flights, we compute and compare 
the statistics of the durations of flights Tf for two trajec- 
tories: a regular quasiperiodic one with the initial posi- 
tion close to the inner border of the vortex core (Fig.[7K) 
and a chaotic one with the initial position close to the 
vortex-core border from the outside (Fig. [7b) . Each full 
rotation of a particle in a frame consists of two flights, 
eastward and westward, with different values of Tf be- 



cause of the zonal asymmetry of the flow. The statis- 
tics for the chaotic trajectory, sticking to the vortex core 
(Fig. [7b), may be considered as a distribution of the dura- 
tions of flights in the vortex-core RIT. The minimal flight 
duration in this RIT is T/ ~ 11 (the flights with smaller 
values of Tf are rare and they occur outside the trap). 
Positions of the local maxima of the PDF for the stick- 
ing trajectory in Fig.[7l3 correlate approximately with the 
corresponding local maxima of the PDF for the regular 
trajectory inside the core in Fig. [TK- The similar corre- 
lations have been found (but not shown here) between 
the local maxima of the PDFs for the lengths of flights 
P{xf) for the interior regular and sticking chaotic tra- 
jectories. These correlations and positions of the peaks 
prove numerically that short flights with \xf\ < 27r and 
11 < Tf < 21 may be caused by the effect of vortex- 
core RIT. We conclude from Fig. [Tb that the vortex- 
core RIT contributes to the statistics of the short flights 
in the range 11 < < 20 for the eastward flights 
with the prominent peak at T/ ~ 11 and in the range 
15 < T/ < 21 for the westward flights with small peaks 
at Tf ~ 17 and 21. 

The effect of the RIT of the secondary-resonance is- 
lands is illustrated in Fig. [H To flnd the characteristic 
times of this RIT we compute two trajectories: a regu- 
lar quasiperiodic one with the initial position inside one 
of these islands and a chaotic one with the initial posi- 
tion close to the outer border of the island. The respec- 
tive PDFs P{Tf), shown in Figs.[8K and b, demonstrate 
strong correlations between the corresponding peaks at 
Tf ~ 12, 23, and 27. Computed (but not shown here) 
PDFs P{xf) for these trajectories conflrm the effect of 
the islands RIT on the statistics of short flights. 

We have shown that the RITs contribute to the short- 
flight statistics, namely, they contribute to all the local 
peaks around Tf ~ 11, 17, and 21 in the short-time PDF 
in Fig. [6l Many other flights, of course, may contribute 
to this statistics. The most abundant among them are 
the flights connected with the hyperbolic saddle trajec- 
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FIG. 8: The secondary-resonance islands trap. The PDFs of 
durations T/ of the eastward (e) and westward (w) flights for 

(a) a regular quasiperiodic trajectory inside the islands and 

(b) a chaotic trajectory sticking to the island's boundary from 
the outside. 

tories and their stable and unstable manifolds. Tracers 
with initial positions close to the stable manifolds are at- 
tracted by the corresponding saddle trajectory perform- 
ing a large number of flights around the points {xs,ys)- 
with the typical values Tf ~ 12 (for the eastward flights) 
and Tf ~ 13 (for the westward flights). The saddle tra- 
jectories, their impact on chaotic transport, evolution of 
their stable and unstable manifolds, and fractal proper- 
ties of chaotic transport will be considered in detail in a 
future publication. 

IV. BALLISTIC-ISLANDS TRAPS 

Besides the rotational islands with particles moving 
around the corresponding elliptic points in the same 
frame, we have found in Paper I ballistic islands situated 
both in the stochastic layer and in the peripheral cur- 
rents. Regular ballistic modes js^ correspond to stable 
quasiperiodic inter-frame motion of particles. Only the 
ballistic islands in the stochastic layer are important for 
chaotic transport. Mapping positions of the regular bal- 
listic trajectories at the moments of time i„ — 2TTn/w 
(n — 1,2,...) onto the flrst frame, we obtain chains 



of ballistic islands both in the northern and southern 
stochastic layers, i. e., between the borders of the north- 
ern (southern) peripheral currents and of the correspond- 
ing vortex cores. A chain with three large ballistic is- 
lands is situated in those stochastic layers. The parti- 
cles, belonging to these islands, move to the west, and 
their mean zonal velocity can be easily calculated to be 
{uf) = -2n/3T = -w/3 ~ -0.0845. There are also 
chains of smaller-size ballistic islands along the very bor- 
der with the peripheral currents. 

We have demonstrated in Paper I a stickiness of chaotic 
trajectories to the borders of those three large baUistic is- 
lands (see Figs. 6 and 7 in Paper I). The Poincare section 
with fragments of two chaotic trajectories in the north- 
ern stochastic layer is shown in Fig.[9K- One particle per- 
forms a long flight sticking to the very border with the 
regular westward current, and another one moves to the 
west sticking to the very boundaries of three large bal- 
listic islands. A magniflcation of a fragment of the bor- 
der and tracks of a sticking trajectory around a smaller- 
size ballistic island are demonstrated in Fig. [9Jd. Fig. 
demonstrates the effective size of the trap of the large bal- 
listic islands with tracks of a sticking trajectory around 
them. 

It is reasonable to suppose that the ballistic-islands 
traps (BIT) contribute, largely, to the statistics of long 
flights with \xf\ >• 2tt. All the baUistic particles, moving 
both to the west and to the east, can flnish a flight and 
make a turn only in the strip shown in Fig. [2l The loci 
of the corresponding turning points have a compHcated 
fractal-like structure that will be studied in detail else- 
where. We consider further only long westward flights, 
taking place in the northern stochastic layer, because it 
is much wider than the stochastic layer between the reg- 
ular central jet and the southern parts of the vortex cores 
where eastward flights take place. 

To distinguish between contributions of the traps of 
different baUistic islands (and, maybe, other zones in the 
phase space) to the statistics of long flights, we com- 
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FIG. 9: (a) Poincare section of the northern stochastic layer 
where stickiness to the very border with the regular westward 
current and to three large ballistic islands are shown. In- 
creased density of points along the border with the peripheral 
current is caused by the traps of the border ballistic islands 
one of which is shown in (b). (c) The trap of the large ballistic 
islands. 

pute for five long chaotic trajectories (up to t = 5 • 10®) 
the distribution of a number of westward flights with 
Tf > 10"^ over the mean zonal velocities (uj) — Xf/Tf 
of the particles performing such flights. The distribu- 
tion in Fig. [To] has a prominent peak centered at the 
mean zonal velocity (u/) ~ —0.0845 which corresponds 
to a large number of long flights of those particles (and 
their trajectories) which stick to the very boundaries 
of the large ballistic islands (see Fig. [9^) moving with 
the mean velocity (u/) ~ —0.0845. The flat left wing 
of the distribution N{{uf)) corresponds to the traps of 
smaller-size ballistic islands nearby the border with the 
peripheral current. There are different families of these 
islands (see one of them in Fig. (SId) with their own val- 




FIG. 10: The distribution of a number of long westward flights 
with Tf > 10^ over their mean zonal velocities (uf). The 
sharp peak corresponds to the trap connected with the very 
boundaries of the large ballistic islands, the left wing — to 
a number of traps of families of the border ballistic islands, 
and the right wing — to the trap situated around the large 
ballistic islands. 

ues of the mean zonal velocity which are in the range 
-0.092 < (uf) < -0.0845. Stickiness to the boundaries 
of the border islands is weaker because they are smaller 
than the large islands and their contribution to the statis- 
tics of long flights is comparatively small. 

The right wing of the distribution N{{uf)) with 
—0.084 < (uf) < —0.075 deserves further investiga- 
tion. The value (w/) ~ —0.075 is a minimal value of 
the zonal velocity for long westward flights possible in 
the northern stochastic layer. Increasing the minimal 
duration of a flight from Tf = 10^ to T/ = (2 ^ 5) • 10^, 
we have found splitting of the broad distribution with 
—0.084 < {uf) < —0.08 into a number of small distinct 
peaks. Comparing trajectories with the values of (u/) 
corresponding to these peaks, we have found that all they 
move around the large ballistic islands. The particles 
with smaller values of (u/) used to penetrate further to 
the south from the islands more frequently than those 
with larger values of (uf) which prefer to spend more 
time in the northern part of the dynamical trap con- 
nected with those islands. Thus, we attribute the right 
wing of the distribution N{{uf)) to an effect of the trap 
situated around the large ballistic islands. 

To estimate the contribution of different BITs to the 
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statistics of long westward flights in Fig. [3] we have com- 
puted the PDFs P{xf) and P{Tf) for particles perform- 
ing westward flights with Xf > 100 and Tf > 1000 
and with the mean zonal velocity (uf) to be chosen 
in three different ranges shown in Fig. [101 —0.092 < 
(uf) < —0.085 (particles sticking to the border islands) 
-0.085 < {uf) < -0.084 (particles sticking to the 
very boundary of three large islands), and —0.084 < 
(""/) S —0.075 (the trap of the three large islands). 
All the PDFs P{xf) decay exponentially but with dif- 
ferent values of the exponents equal to ~ —0.005 and 
1/ ~ -0.0018 ~ -0.0014 for the traps of border and the 
large baUistic islands, respectively. The tail of the PDF 
P{xf) for westward flights, shown in Fig. [H decays ex- 
ponentially with ly ~ —0.0014. Thus, the contribution 
of the large island's BIT to the statistics of long west- 
ward flights is dominant. As to temporal PDFs P{Tf) 
for westward long flights, they are neither exponential 
nor power-law like with strong oscillations at the very 
tails. The slope for the border BITs is again smaller 
than for the large ballistic islands trap. 

V. CONCLUSION 

A meandering jet is a fundamental structure in oceanic 
and atmospheric flows. We described statistical proper- 
ties of chaotic mixing and transport of passive particles in 
a kinematic model of a meandering jet flow in terms of dy- 
namical traps in the phase (physical) space. The bound- 
aries of rotational islands (including the vortex cores) in 
circulation zones are dynamical traps (RITs) contribut- 
ing mainly to the statistics of short flights with 1 < 27r. 
Characteristic times and spatial scales of the RITs have 
been shown to correlate with the PDFs for the lengths 
Xf and durations Tf of short flights. The boundaries of 
ballistic islands in the stochastic layers (including those 
situated along the border with the peripheral current) 
are dynamical traps (BITs) contributing mainly to the 
statistics of very long flights with ^ 27r. 



Dynamical traps are robust structures in the phase 
space of dynamical systems in the sense that they present 
at practically all values of the corresponding control pa- 
rameters. We never know exact values of the parameters 
in real flows, especially, in geophysical ones. We don't 
know exactly the structure of the corresponding phase 
space, however, we know that typical features, like is- 
lands of regular motion, vortices, and jets, exist in real 
flows (see their images in some laboratory flows in Ref. 
[sii). In this paper we chose speciflc values of the control 
parameters for which speciflc PDFs have been computed 
and explained by the effect of those dynamical traps that 
exist under the chosen parameters. We have carried out 
computer experiments with different values of the control 
parameters and found that the phase space structure has 
been changed, of course, with changing the values of the 
parameters, but the corresponding RITs and BITs with 
speciflc temporal and spatial characteristics have been 
found to contribute to the corresponding statistics. 

There are, of course, another structures in the phase 
space that may have an effect on chaotic transport of 
passive particles. The most important of them are hy- 
perbolic saddle trajectories and their stable and unsta- 
ble manifolds. They should contribute to the statistics 
of transport on all the scales. We leave elaboration on 
this topic and on the topics of fractality of chaotic mix- 
ing and transport in the meandering jet flow to a future 
publication. 

After flnishing our work, we were acquainted with Ref. 
[3^ where meridional chaotic transport, associated with 
a similar kinematic model of a meandering jet, has been 
studied by the method of lobe dynamics 33|- In differ- 
ence from our study of zonal chaotic transport, a geomet- 
ric structure of cross jet transport has been considered 
in Ref. where values of the control parameters have 
been chosen to be sufficiently large to break up the cen- 
tral jet as a barrier to transport of particles across the 
jet. The mechanisms for particles to cross the jet have 
been described in terms of lobe dynamics and the mean 
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time to cross the jet for particles entering the jet and the 
mean residence time for particles in the jet have been 
estimated in Ref. [34|. We have studied a more realistic 
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situation (at least, for surface oceanic jet currents) when 
the jet is an absolute barrier to cross jet transport and we 
explained statistical properties of transport in terms of 
dynamical traps of rotational and ballistic islands. The 
method of lobe dynamics is hardly applicable for study 
zonal chaotic transport since it is practically impossible 
to trace out lobe evolution for a large number of frames. 



The work was supported by the Russian Foundation for 
Basic Research (Grant no. 06-05-96032), by the Program 
"Mathematical Methods in Nonlinear Dynamics" of the 
Russian Academy of Sciences and by the Program for 
Basic Research of the Far Eastern Division of the Russian 
Academy of Sciences. 



[1] S.V. Plants, M.V. Budyansky, M.Yu. Uleysky, and 

G.M. Zaslavsky, Chaos 16, 033117 (2006). [19 
[2] G.M. Zaslavsky, Phys. D. 168-169, 292 (2002). 
[3] H.A. Dijkstra, Nonlinear physical oceanography (Dor- [20 

drecht, Kluwer, 2000). 
[4] S. Wiggins, Annu. Rev. Fluid Mech. 37, 295 (2005). [21 
[5] R.T. Pierrehumbert, Phys. Fluids 3, 1250 (1991). 
6] M. Cencini, G. Lacorata, A. Vulpiani, and E. Zam- [22 
bianchi, J. Phys. Oceanogr. 29, 2578 (1999). [23 
7] T.F. Shuckburgh and P.H. Haynes, Phys. Fluids 15, 3342 

(2003). [24 
8] V.I. Arnold, C. R. Hebd. Seances Acad. Sci. 261, 17 [25 

(1965) . 

9] M. Henon, C. R. Hebd. Seances Acad. Sci. 262, 312 [26 

(1966) . 

0] H. Lamb, Hydrodynamics (Dover, New York, 1945). [27 
1] G.M. Zaslavsky, Phys. Rep. 371, 461 (2002). [28 
2] G.M. Zaslavsky, Hamiltonian chaos and fractional dy- 
namics (Oxford University Press, Oxford, 2005). [29 
3] D. del-Castillo-Negrete and P.J. Morrison, Phys. Fluids 

A 5, 948 (1993). [30 
4] R.M. Samelson, J. Phys. Oceanogr. 22, 431 (1992). [31 
5] K. Ngan and T. Shepherd, J. Fluid Mech. 334, 315 
(1997). 

6] G.C. Yuan, L.J. Pratt, and C.K.R.T. Jones, Dyn. Atmos. [32 
Oceans. 35, 41 (2002). [33 
7] M.V. Budyansky, M.Yu. Uleysky, and S.V. Prants, Ne- 

lineinaya Dinamika 2, 165 (2006) [in Russian]. 
8] J. Sommeria, S.D. Meyers, and H.L. Swinney, Nature 



(London) 337, 58 (1989). 

R.P Behringer, S.D. Meyers, and H.L. Swinney, Phys. 
Fluids A 3, 1243 (1991). 

T.H. Solomon, E.R. Weeks, and H.L. Swinney, Phys. 
Rev. Lett. 71, 3975 (1993). 

T.H. Solomon, E.R. Weeks, and H.L. Swinney, Phys. D 
76, 70 (1994). 

A.S. Bower, J. Phys. Oceanogr. 21, 173 (1989). 

A. S. Bower andH.T. Rossby, J. Phys. Oceanogr. 19, 1177 
(1989). 

B. V. Chirikov, Phys. Rep. 52, 263 (1979). 

R.S. Mackay, J.D. Meiss, and I.C. Percival, Phys. D. 13, 
55 (1984). 

K.V. Koshel and S.V. Prants, Physics-Uspekhi 49, 1151 

(2006) [Uspekhi Fizicheskikh Nauk 176, 1177 (2006)]. 

G.M. Zaslavsky, Chaos 5, 653 (1995). 

G.M. Zaslavsky and B. Niyazov, Phys. Rep. 283, 73 

(1997). 

B.A. Petrovichev, A.V. Rogalsky, R.Z. Sagdeev, and 
G.M. Zaslavsky, Phys. Lett. A 150, 391 (1990). 
V. Rom-Kedar and G. Zaslavsky, Chaos 9, 697 (1999). 
J.M. Ottino, The kinematics of mixing: stretching, 
chaos, and transport (Cambridge University Press, Cam- 
bridge, 1989). 

F. Raynal and S. Wiggins, Phys. D 223, 7 (2006). 
W.S. Wiggins, Chaotic Transport in Dynamical System 
(Springer- Verlag, New York, 1992). 



